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Abstract 

A recently proposed molecular model is discussed as a non-trivial extension of 
the Ising model. For d = 2 the two models are shown to be equivalent, while 
for d > 2 the molecular model describes a peculiar second order transition 
from an isotropic high temperature phase to a low-dimensional anisotropic 
low temperature state. The general mean field analysis is compared with the 
results achieved by a variational Migdal-Kadanoff real space renormalization 
group method and by standard Monte Carlo sampling for d = 3. By finite 
size scaling the critical exponent has been found to be = 0.44 it 0.02 thus 
establishing that the molecular model does not belong to the universality class 
of the Ising model for d > 2. 

PACS numbers: 64.60.Cn, 64.60.Fr, 62.50.-hp, 05.50.-Fq 
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I. INTRODUCTION 



The molecular model has been first proposed |l| as a very simple (i-dimensional lattice 
model which incorporates some degree of frustration and thus describes some aspects of 
molecular orientation in covalently bound molecular solids. Most molecular liquids retain 
their molecular structure even in the solid phase, where some long range order usually shows 
up as a consequence of inter-molecular interaction. However in the solid the orientational 
order of the molecules may change according to the thermodynamic conditions giving rise 
to quite rich phase diagrams as recently observed for hydrogen under high pressure [Q. 
Moreover orientational ordering is responsible for several phase transitions occurring even 
in the liquid phase (liquid crystals), and lattice models have been found to be useful for 
describing such transitions. 

The most studied models describe a molecular interaction which arises from dipole fluc- 
tuation, is weak and gives rise to the observed three-dimensional ordering of most molecular 
Van der Waals solids. The thermodynamic behaviour of such weakly interacting systems 
can be analyzed in terms of 0(3) symmetric vectorial models. 

Conversely the molecular model was motivated by a description of the almost-covalent 
molecular solids where the interaction has a covalent main component and is characterized 
by some level of frustration (since the coordination number for the covalent bond is quite 
low). In such solids each molecule must choose a few partners and cannot accept any further 
invitation. The lower is the allowed coordination number, the higher the frustration which 
gives rise to the low- dimensional structures observed in polymers (one-dimensional) or in 
iodine and hydrogen halides (two-dimensional). In iodine, where the covalent nature 
of the interaction is out of doubt even one-dimensional zigzag chains have been reported 
under pressure Q. Moreover we expect that a covalent interaction should show up for all 
the molecular solids under high pressure, as the inter-molecular distance approaches the 
intra-molecular length, provided that some important structural transition does not occur 
first (like dissociation). Quite recently the existence of charge transfer between O2 molecules 
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has also been reported under pressure 

The molecular model is a simple frustrated lattice model which can describe some aspects 
of a covalently bound molecular solid. It consists of a d-dimensional hypercubic lattice with 
a randomly oriented linear molecule at each site. In its simplest version each molecule is 
only allowed to be oriented towards one of its nearest neighbours. There is an energy gain 
for any pair of neighbours which are oriented along their common bond (a covalent bond). 
The existence of preferred orientational axes breaks the rotational invariance of the single 
molecule as it is likely to occur for any real molecular system under pressure. In fact even 



in hydrogen, the broken symmetry phase transition which is observed under pressure |10 



has been recently shown |jTT| to be affected by the presence of a crystal field which breaks 
the isotropy. 

Quite recently, similar lattice models have been used for describing the diffusion of par- 
ticles and molecules inside a polymer, and the growth of one-dimensional islands (polymeric 
chains) [jl2[. The molecular model has already stimulated some recent work on molecular 
orientation in nitrogen |TB|JT^ which goes back to the phenomenology laid down by Pople 



and Karasz [|rH|. It had been argued [|T7| that the weak intermolecular bonds between two 



N2 molecules should be favourable for the formation of an orientationally disordered "plas- 
tic crystal" solid phase, and should lead to freezing into an orientationally ordered phase. 



More recent experimental data on nitrogen [16| confirm the existence of an orientational 



disordering temperature in the solid below the melting temperature. However, as far as we 
know, the molecular systems which are more closely described by the molecular model are 
the hydrogen halides HX{X = F, CI, Br, I). Their low-temperature structures are known 
to consist of planar chains of molecules in the condensed state while a totally disordered 
structure has been observed with increasing temperature at ambient pressure 0] . Moreover 
the opposite transition, from orientational disorder to an ordered chain structure, has also 
been reported by increasing pressure 0. 

The molecular model undergoes a transition from an high-temperature (or weakly in- 
teracting) fully isotropic disordered system, to a low-temperature (or strongly interacting) 



anisotropic low dimensional broken-symmetry phase. As a consequence of frustration the 
breaking of symmetry is accompanied by a sort of decomposition of the system in low- 
dimensional almost independent parts, as observed in solid iodine and hydrogen halides. 
Such remarkable behaviour requires a space dimension d > 2, while for d = 2 the model is 
shown to be equivalent to the exactly solvable two-dimensional Ising model [|T^. As shown 
by Monte Carlo calculations, in the broken-symmetry phase the system displays the presence 
of correlated chains of molecules (polymers) which point towards a common direction inside 
each two-dimensional sub-set of the lattice (plane). Such planes are weakly correlated in the 
low-temperature phase, and the system has a two-dimensional behaviour even for d > 3. 

In this paper the relevance of the molecular model as a non-trivial extension of the 
Ising model is pointed out. Thus, apart from the physical motivations, the model is fully 
examined and the phase transition is described by several methods: mean-field, real space 
renormalization group and numerical simulations. Exactly solvable models are important for 
our understanding of more complex systems, and provide a test for approximate techniques. 
The (i- dimensional molecular model shares with the Ising model the d = 2 realization, since 
their equivalence for d = 2 has been proven to be exact |I| . In this paper we will focus on the 
d = 3 model, but we will take advantage of the existence of an exactly solvable realization 
for d = 2. For > 2, as the frustration increases, the model shows a very different behaviour 



compared to the Ising or Potts models. These last show a fully (i- dimensional broken- 
symmetry phase while the molecular model is characterized by a low-dimensional ordering 
inside the planes with negligible correlation among different planes. Moreover for d = 3 
the molecular model is shown to belong to a different universality class, since its critical 
exponent v turns out iohe v = 0.44 ± 0.02 by finite size scaling. We expect that such new 
universality class should describe a broad group of isotropic physical systems characterized 
by a low- dimensional ordering in their low-temperature phase. Such broad class of phase 
transitions should be explored by experiments in order to compare with the theoretical 
predictions for the critical universal properties. In such respect the driving parameter does 
not need to be the temperature, as the bond strength can be directly modified by a change 



of pressure in several systems. 

This paper is organized as follows. Section II contains a formal definition of the d- 
dimensional molecular model, and a proof of its equivalence to the Ising model for d = 2. 
In Section III the mean-field solution is discussed for the generic dimensional model. In 
Section IV a modified variational Migdal-Kadanoff method is presented and its application 
is discussed for d = 2 and d = 3. At variance with a previous calculation I^Ol which yielded 
a quite poor result, the variational method is shown to work very well provided that some 
assumptions are made on the nature of the broken phase. Section V contains the results 
of a numerical simulation by Monte Carlo sampling, and the numerical estimation of both 
critical temperature and exponent by finite size scaling. In Section VI the main findings are 
summarized and discussed. 

II. THE MOLECULAR MODEL 

The molecular model was first described in Ref. |Q. We briefly describe its formal 
definition in order to fix the notation. Let us consider a dimensional hypercubic lattice, 
with a randomly oriented linear molecule at each site. The molecules are supposed to be 
symmetric with respect to their centre of mass which is fixed at the lattice site. Only a 
discrete number of space orientations are allowed for each molecule: we assume that each 
of them must point towards one of its 2d first neighbour sites. This choice can be justified 
by the existence of covalent interactions along preferred axes. Then each molecule has d 
different states corresponding to molecular orientation along the hypercube axes (molecules 
are symmetric). Finally, each couple of first neighbour molecules, when pointing the one 
towards the other, are assumed to gain a bonding energy for their directional covalent bond 
(they touch each other). As shown in Fig.l for d = 2, bonding in a direction excludes any 
possible bond along the other {d— 1) directions. The coordination number is 2 for any value 
of d, and the frustration increases with increasing d. 

According to such description we introduce a versor variable Wj. for each of the sites 
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r of the lattice, with Wr G {xi, £2, • ■ • Xd} pointing towards one of the d hypercube axes Xa- 
The versors x^ are assumed to be orthonormal: x^-x^ — S^-y- The partition function follows 

Z = ^e^ = Eexp 

{w} {w} 

where {tv} indicates a sum over all configurations, a runs from 1 to d, and the lattice spacing 
is set to unity. The inverse temperature (3 (in units of binding energy) can be negative for 
a repulsive model, but is assumed positive in the molecular context. 

The model may be generalized by introducing an external rf-dimensional vectorial field 
h{a) at each link. The dependence on a means that the field differs according to the space 
direction a of the lattice link which joins the sites. The modified partition function reads 

Zh^^e^'^ =J2 CXP ] 4/3 J2 [(^r • Xa) (Wr+xo. ' ^a) + h(Q;) • Wrr + ^{o) ' t&r+x^] \ (2) 
{w} {w} I r,a J 

It is evident that if the field satisfies the condition 

Eh(«) = (3) 

a 

then Sh does not depend on h and Sh = -S". In such case the extra degree of freedom provided 
by h can be regarded as a sort of internal symmetry of the model. This global symmetry 
can be made local by allowing the field h to depend on the site position r. We will only take 
advantage of the global symmetry in this paper. We notice that such symmetry cannot be 
seen as a gauge invariance, since in lattice gauge models any gauge change leaves the energy 
gain unchanged at any link. Here the field h changes the energy gain of all the links while 
the whole action is invariant. 

Adopting a more compact notation, the partition function reads 

Zn = Y.^'^ = Y.^^''-''^^''''^ (4) 
{&} {&} 

where the Lagrangian density C follows as 

£(r, a) = wlMa {f3, h)wr+x„ (5) 



4/5 V(Wr • Xa){Wr+x^ " £«) 



(1) 
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Here the canonical d-dimensional column vector representation of is employed with xi = 
(1, 0, . . .), £2 = (0, 1, . . .), etc. The dxd matrix does not depend on the configurations 
of the system, and entirely characterizes the model. 

The global symmetry of the action provides a simple way to show the equivalence between 
molecular and Ising models for d = 2. For the two-dimensional lattice the condition is 
satisfied by the field h(l) = h{xi — X2), h(2) = — h(l). The matrix Ma follows 

/4(3{l + 2h) 

Mi = 

V -8(3h^ 
( -8(3h \ 

V 4/?(i + 2/i)y 

Then for h = -1/4, Mi = M2, and C reads 

L(r, a) = (3 + wl\ Wr+x„. (7) 

Identifying the two-dimensional column versors w with spin variables, apart from an inessen- 
tial factor, Z reduces to the partition function of a two-dimensional Ising model 



2I3N 



■ Zising (8) 



and is exactly solvable. For (3 +00 a ground state is approached with all the molecules 
oriented along the same direction, and with formation of one-dimensional polymeric chains 
(Fig. 2a); for P —00 the repulsive model approaches a zero-energy (no bonds) ground 
state analogous to the antiferromagnetic configuration of the Ising model (Fig.2b). 

For d > 3 the analogy with the Ising model breaks down, and this is evident from a 
simple analysis of the ground state configuration. Due to frustration the model has an 
infinitely degenerate ground state in the thermodynamic limit N 00. For instance, in the 
case d = 3, the minimum energy is obtained by orienting all the molecules along a common 
direction, as for d = 2. However the ground state configuration is not unique: the number 
of molecular bonds does not change if we rotate together all the molecules belonging to 
an entire layer which is parallel to the original direction of orientation. As a consequence 



of frustration the total degeneration is 3(2^^^ ^)), and the system could even behave like 
a glass for the large energy barriers which separate each minimum from the other. The 
phase diagram is expected to be quite rich, with at least a transition point between the high 
temperature disordered phase and an ordered broken-symmetry low temperature phase. 



III. MEAN FIELD APPROXIMATION 

For the generic (i-dimensional model, some analytical results can be obtained in Mean- 
Field (MF) approximation: neglecting second order fluctuation terms 

where Aq, = {wj. ■ Xa) is an average over the configurations, and J2a = 1 (with the 
obvious bounds < Aq < 1). Here the order parameter A^ gives the probability of finding 
a molecule oriented along the direction of Xa- The partition function factorizes as 

N 



(10) 



and the free energy follows 



FmF = -^^OgZMF = 4E - ^log (E^'"^'^) • (11) 



The derivative with respect to A„ yields, for the stationary points 



A. = ^^-^sm:: (12) 



which satisfies the condition J2a = 1- 

In the high temperature limit (3^0 eq.(|T2|) has the unique solution A^ = 1/d which 
reflects the complete random orientation of molecules. In the opposite limit /3 — oo, apart 
from such solution, eq. ([T^) is satisfied by the broken-symmetry field A^ = 1, A^ = 
for a ^ fi, which obviously corresponds to a minimum for Fmf- Then at a critical point 
P = Pc the high temperature solution must become unstable towards a multivalued minimum 
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configuration. The Hessian matrix is easily evaluated at the stationary points by using 
eqs.(0) and (|n|): 

= = ~ ^^^'^ + ^^^^'^^ ^^^^ 

In the high temperature phase (/3 < /3c), inserting = ^/d, the eigenvalue problem 

det \H^, - \S^,\ = (14) 

yields 

(l- f-^)' '■(1-A) = 0. (15) 

Thus the Hessian matrix is positive defined if and only if A = {1 — 8(3/d) > 0. Beyond the 
critical point (3 = (3c = {d/8) the solution = l/rf is not a minimum, and a multivalued 
minimum configuration shows up. Such result obviously agrees with the MF prediction for 
the Ising model, Rising = l/(2rf), only for the special dimension d = 2. For > 2 we observe 
an increase of f3c with d, to be compared to the opposite trend shown by the Ising model. 
Such behaviour may be interpreted in terms of the low dimensionality of the ordered phase. 
Due to frustration the ordering may only occur on a low dimensional scale: for instance in 
three dimensions each layer has an independent internal ordering. Thus we expect a larger 
/3c for d > 2 since the increasing of d only introduces larger fiuctuations, with each molecule 
having {d — 2) allowed out-of-plane orientations. For d = 3 the low temperature phase can 
be regarded as a quenched disordered superposition of layers which are internally ordered 
along different in-plane directions. As a consequence of frustration the system shows a two- 
dimensional character below the critical point while behaving as truly three-dimensional in 
the high temperature domain. In MF the neglecting of some fiuctuations usually leads to a 
critical temperature which overestimate the exact value (i.e. the critical inverse temperature 
(3c is underestimated). For d = 2 the MF prediction is (3c = 0.25 to be compared with the 
exact value (3c = 0.4407. For d = 3 the MF prediction (3c = d/8 = 0.375 should provide a 
lower bound to the unknown exact value. 
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IV. VARIATIONAL MIGDAL-KADANOFF APPROXIMATION 



According to the Migdal-Kadanoff [18,21| method, a hnk displacement may be introduced 
by considering that the configurational average of the Lagrangian density L in eq.(|) must 
be translationally invariant 

(/:(r,a)) = (£(r',a)) (16) 

then defining 

r„(r,r') = /:(r,a)-/:(r',a) (17) 

we can state that (rQ,(r, r')) = and the same holds for any sum F over an arbitrary set of 
such terms 

r = 5:r.(r,rO. (18) 

Replacing the action Sh by the sum Sh + T, and assuming that the condition (|^) is verified 
(so that we can drop the h in S^, and which are invariant), the modified partition function 
Zp can be approximated by cumulant expansion as 

Zr = E e^+r = Z ■ I/) ^ Z [e<r> ■ eK(r^)-<r>^)l ^^g^ 
{^} ^ ^ 

then, since (F) = 0, 

Zr^Z-e5<r'> (20) 



For instance, the sum in equation (|T8|) could run over all a 7^ 1, and for appropriate values 
of the vectors r, r', in order to yield a displacement of links which are orthogonal to X\. To 
second order in F, the error introduced by link displacement is controlled by the exponential 
factor in equation (pll). 

Link displacement breaks the internal symmetry of the model, so that Z^ is no longer 
invariant for any field change subject to the condition (^. Then we may improve the 
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approximation by using the extra freedom on the choice of h for minimizing the difference 
between the approximate partition function Zy and the exact Z. 

If h satisfies the condition (|^) then ah satisfies such condition as well for any choice of 
the scalar parameter a. Then a special class of invariance transformations can be described 
by a change of the strength parameter h, assuming the field h as proportional to h. The 
following discussion could be easily generalized to other classes of transformations described 
by more than one parameter. Since F is linear in the field h, then in general 

= [A + hBf (21) 

where A and B depend on the configuration of the system. For the average we have 

(r^) = {A^) + 2h{AB) + h''{B^) (22) 

This last equation, inserted in eq. (pOD leads to the following considerations: i) the coefficient 
(-B^) is positive defined, then the average (F^) always has a minimum for an appropriate 
value of h = ho] ii) in general (AB) ^ then ^ 0, and a direct use of the Migdal- 
Kadanoff method on the original model (with no field considered) would yield a larger error; 
iii) to the considered order of approximation Z^ is stationary at h = h^, and is symmetric 
around that point, then all the physical properties described by such partition function must 
result symmetric with respect to Hq. Moreover, at the same order of approximation, any 
physical observable / will acquire an unphysical dependence on h, and the symmetry around 
ho requires that ^ = for /i = /iq- Then we expect that all such observables should be 
stationary at h = h^, and their best estimate should coincide with the extreme value. 

As a consequence of the above statements, the Migdal-Kadanoff method can be improved 
by taking advantage of the global symmetry of the model. By use of the approximate 
partition function Z^ the critical temperature acquires a non-physical field dependence, but 
the best estimate of jSc is its stationary value corresponding to h = h^. The method can be 
seen as a variational method with the best approximation achieved by the minimum in the 
inverse temperature. 
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Such stationary condition resembles the principle of "minimum sensitivity" introduced 
by Stevenson for determining the best renormalization parameters whenever the phys- 
ical amplitudes depend on them (and they should not). In our context, since the critical 
temperature should not depend on the choice of the field strength h, the best value for 
such field is the one which makes the critical temperature less sensitive i.e. the stationary 
point. However, according to equation (|20|) and (p2|), here we have got a formal proof of the 



stationary condition up to second order of the cumulant expansion. 

The method may be used by performing a displacement of links that are orthogonal to 
Xi, and then a one-dimensional decimation along the a = 1 axis. According to such program 
let us define the alternative d x d matrix ta{l3, h) 

e^^^'"^ = w%{P,h)w,+^^ (23) 

The partition function follows 

Zh = Y.X{W^(P^h)w,+-<] (24) 

{w} r,a 

After link displacement and decimation along the a = 1 axis, the modified partition function 
reads 

ZT = Y.X{[^lU^M^r+,J\ (25) 

{w} r,a 

where the sum and the product run over the configurations and the sites of the new decimated 
lattice, and 

ti(/3,/i) = [ti(/3,/i)]' (26) 

t„(/5, h) = taiXp, h) for a ^ 1 (27) 

with A being the scale factor between the new and the old lattice. A renormalized inverse 
temperature [3a may be defined according to 

ti(A/i) = ti(/3i,/i) (28) 
12 



= \(3 for a ^ 1 (29) 

Eventually, the same scaling operation should be performed consecutively for all the di- 
rections in order to obtain an hyper-cubic lattice again. For any finite scaling parameter 
A > 1 the renormalized inverse temperature is anisotropic, but an isotropic fixed-point can 
be recovered in the limit A ^ 1. The equations (^8D , (|29|) define the flow of the renormalized 
inverse temperature, which changes for any different value of the field strength h. Equation 
(p8|) has a more explicit aspect in the representation of the common eigenvectors of the 
matrices ti and ti = [ti]^. The rank of such matrices is 2 for any space dimension d, as can 
be expected from the definition of the model. Then both the matrices can be represented in 
terms of the two no n- vanishing eigenvalues rji, ri2, which are functions of (3 and h. Assuming 
that 7^ 0, and defining 

f{P,h) = ^ (30) 

apart from a regular multiplicative factor for the partition function, the scaling equation 
(^) reads 



[f{(3,h)]' = f0„h) (31) 

For any fixed h, the fixed points follow through the standard Migdal-Kadanoff equations 

fiX^-'(3^,h)]' = fiX^-'(3^,h). (32) 

When A is analytically continued up to 1 such equations give the same isotropic fixed point 
Pc- In fact, the expansion of equations ( p2D around A = 1 implies (up to first order in A — 1) 



\nf{P,,h) = -{d-l)P, 



ld[_ 
fdp 



(33) 

/3c 



which is an implicit equation for /3c- Such equations yield their best estimate of jSc when the 
strength of the field h is set to the stationary value Hq. 

It is instructive to evaluate the stationary point Hq for the case d = 2 which is equivalent 
to the two-dimensional Ising model for the choice h = hj = —1/4, as shown in section II. The 
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h invariance of the exact partition function guarantees the equivalence of the two models 
for any choice of h hj. However, the mere application of the Migdal-Kadanoff equations 



to predict the exact fixed point in the limit A 1. That can also be checked by inserting in 
equation ( ^3|) the exact expression for the fixed point of the two-dimensional Ising model. 
Such contradictory results are not surprising since, as already discussed, link displacement 
breaks the h invariance of the model, and the approximate solution is thus dependent on 
the choice of h. We would like to test the variational method on this exactly solvable model: 
we look for the stationary point of the function f{(3,h). The matrices follow from the 
equations (H) 



It can be easily shown that if the derivative of / is zero at a given h independent of /3, then 
the solution jSc of Eq.(p3|) is stationary at that h value. Differentiating with respect to x, 
we find that the derivative of / vanishes for x^ = 1/6, which yields h = —1/4 = hj for 
any jS. As expected, this is the required value in order to recover the Ising model. Thus 
the Migdal-Kadanoff approximation gives an improving estimate of the critical point as we 
move from the molecular towards the Ising representation (where the approximation yields 
the exact fixed point). We stress that all such representations are equivalent due to the h 
invariance of the action. 

For > 2 no equivalence to standard studied models has been found, and the behaviour 
seems to be dictated by the strong frustration which does not allow an higher coordination 



(p^) to the simple h = molecular model fails to predict even the existence of the fixed 
point. On the other hand, for h = hj, the very same recurrence equations (^) are known 




(34) 



where b = exp{Ai3), and x = exp{Aj3h). Then for the eigenvalues we obtain 




(35) 
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number than two, even for higher dimensions. We will focus on the d = 3 model in order to 
compare the results with the Monte Carlo findings of the next section. First of all the fields 
h{a) must be defined. An isotropic choice would be 

h(l) = h{xi - - ^xs) 



h(2) = h{x2 - - ^xi) 
h(3) = h{x3 - ixi - ix2) 



(36) 



The matrix ti follows 





/ g4/3+8/3/i 


g2/3h 




/ bx'^ 


\fx 




h = 






g-4/3/i ^ 


= y/x 


l/x 


l/x 






g-4/3h 


g-4/3/i ^ 




l/x 


l/x) 



(37) 



Then from the eigenvalues we obtain 



/(/5,/^) = ^ 



(bx^ + 2) - ^{bx^ - 2)2 + 8x3 



(38) 



V2 {bx^ + 2) + ^(6x3 - 2)2 + 8x3 

Differentiating with respect to x we find that the derivative of / only vanishes for x^ = 2/b. 
This is equivalent to say 

1 1 



In 2 



(39) 



3 up 

which depends on /?. For such field strength hm{P) the ratio between the eigenvalues reduces 



to 



/(/?,/iJ=tanh(/3) 



(40) 



which is exactly the same expression holding for the Ising model |T^. However we must 
point out that in such case hm is not the stationary point Hq. Since depends on P, the 
vanishing of the derivative of / does not imply that the solution of eq.(|33|) is stationary. In 
fact, for d = 3, the choice h = hm yields the known poor result I3c = 0.1398 by insertion 
of Eq.(^Dp in Eq. (P^D . On the other hand, by insertion of the general expression for / 
Eg . (^8l) , the scaling equation ( p3|) can be numerically solved for f3c as a function of h. At the 
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stationary point (3c has a minimum, and thus the variational method yields an even worse 
prediction (P^ ~ 0.12 at the stationary point). These shortcomings show that the isotropic 
d = 3 variational method does not suite the molecular model. Actually both MF and Monte 
Carlo methods predict a larger /3c and, as pointed out at the end of the previous section, 
the exact Pc should be larger than the MF prediction Pmf — 0.375. 

We could have guessed such disagreement since we are using an isotropic version of the 
variational Migdal-Kadanoff method for a system which is not isotropic in its ordered phase. 
At the transition point the system choices a direction, as is usual for any symmetry breaking 
mechanism. However, at variance with usual models, in the ordered phase the correlation 
length cannot be isotropic: order occurs inside all layers which are orthogonal to the chosen 
direction, while there is a negligible correlation along such direction. It would be more 
sensible to describe the ordering which takes place inside a single layer, thus neglecting 
any correlation among different layers. Inside each layer the correlation length is isotropic, 
and the d — 2 variational Migdal-Kadanoff method should give a better description of the 
transition. The same argument should hold for the generic d-dimensional molecular model. 
Moreover, despite the cost of this further approximation, the Migdal-Kadanoff method is 
known to work better for the lower dimensions, and ad — 2 variational method could provide 
a tool for describing the generic d-dimensional molecular model even for d > 3. 

A d — 2 version of the variational method requires a different choice for the fields h{a) 
which do not need to be isotropic any more. Let us take the same field we used in section 
II, namely h(l) = h{xi — £2), h(2) = — h(l) and h(3) = 0. The matrix ti follows 



1 1/x 
V X 1/x 1 J 



(41) 



U = 1 e-^^'^ e-^^'^ 
V e^^'^ e-^^'^ 1 / 
Notice that this is a 3 x 3 matrix since we are using the d — 2 method but we are still dealing 
with ad — 3 molecular model. The two matrices ti and t2 share the same eigenvalues. Their 
ratio is 
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(bx^ + 1 + - J(bx^ - 1 - l/x2)2 + 4(1 + X2) 

h) = ^ ^ ^ V / ^ ^ ^ (42) 

(6x2 + 1 + l/2;2) + y^(6x2-l-l/x2)2+4(l + x2) 

Inserting this result in the scahng equation (|33|) evaluated at d = 2 yields an implicit 
equation for /3c versus /i. The numerical solutions are reported in Fig. 3. They share most 
of the features of the d = 2 molecular model: (i) There are several solutions but there is 
no repulsive fixed point for h = 0; (ii) The physical solution starts at a negative h which in 
this case is h ^ —0.226; (iii) The physical solution has just one stationary point Hq where 
(3c reaches its minimum value. However in this case the stationary point is at ho = —0.2349 
where f3c = 0.6122. This best estimate of the critical point is not too far from the finite 
size scaling prediction of the next section jSc = 0.53. The result is encouraging, and gives 
us more confidence in our understanding of the physics described by the molecular model. 
Strictly speaking, this d = 2 variational method describes the transition occurring in a single 
layer of molecules. However, at variance with the d = 2 molecular model, each molecule 
is now allowed to be oriented along three different axes (two in-plane and one out-of-plane 
orientations). Thus this reasonable prediction for jSc could be regarded as an indirect proof 
that the correlation between two different layers is negligible, and that in the ordered phase 
the system behaves as a truly two-dimensional one. 



V. MONTE CARLO SAMPLING 

In order to check the prediction achieved by different approximate methods it would be 
desirable to have an accurate numerical estimate of the critical temperature. That can be 
easily obtained by finite size scaling. Moreover, according to the scaling hypothesis, the 
critical exponent u can be extracted by the numerical data with a good accuracy. 

Cubic samples N x N x N with = 10, 15, 20, 25, 30 have been considered. All the 
averages have been evaluated by Monte Carlo sampling with no special boundary conditions. 

In this model any ordering is characterized by the presence of some degree of correlation 
along one-dimensional chains of molecules. For d = 3 there are 3 x N x N different chains 
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in each sample. Each chain may be labelled by its direction a — 1,2,3 and by a couple 
of integer coordinates /i, I2 running over a lattice layer orthogonal to the axis x^- For any 
chain we define an order parameter 

1 ^ 

m{a, h, h) = ^ b^iJoc, h, h) ■ (43) 

where w{Ja,Ii,l2) is the versor Wr at the chain site r whose integer coordinates are de- 
termined by Ja along the chain and by the couple /i,/2 in the orthogonal directions. If 
there is no correlation at all (/3 — > 0) then m{a, 7i, I2) ~ 1/3 for any chain in the sample. 
By averaging over all the chains of each sample and over all the configurations, we obtain 
(m) = 1/3. For large N, according to the central limit theorem, in this statistical ensemble 
the variable m follows a gaussian distribution centered at its average value. In the opposite 
limit {f3 — > 00) a third of the chains in each sample have a large m ~ 1, while m ~ for 
two thirds of them. Since any intermediate value of m is unlikely, the statistical distribution 
of m can be regarded as the superposition of two different peaked distributions centered at 
m = and m = 1. If is large enough, and for a large number of configurations, such 
distributions are very peaked and their width is very small. Actually, just below the critical 
point the gaussian distribution already splits in a double-peak distribution. We can monitor 
the transition by use of the new variable 7 

7 = / = (44) 

J {w?) — (m)^ 

By its definition the configurational average of 7 is vanishing (7) = and the second moment 
(7^) = 1. The variable 7 only differs from m for a shift and a rescaling, thus the statistical 
distribution for 7 follows the same trend already discussed for m. However the fourth 
moment s = (7"^) is now strongly dependent on the number of peaks characterizing the 
statistical distribution. For a single gaussian s = 3 exactly. Below the critical temperature 
the distribution splits. In the thermodynamic limit N ^ 00 the width of each peak vanishes, 
while the two peaks separate by a finite quantity. For instance assume that just below the 
critical point a third of the chains yield m ~ 1/3 -|- e where e is a very small increase in the 
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chain correlation which breaks the symmetry of the sample. The other two thirds of chains 
must yield mRil/3 — e/2 since by its definition (m) = 1/3 exactly. Neglecting the width 
of the peaks we may approximate the statistical distribution for m as the superposition of 
two delta-functions with weight factors: 

P{m) = ^5(m - ^ + |) + - ^ - (45) 

By use of such approximate statistical distribution the calculation of the fourth moment 
s = (7"^) is straightforward and gives s — 1.5 for any e, no matter how small. This is 
one half of the single gaussian value. Thus in the thermodynamic limit we expect that the 
fourth moment s should behave like a step function with constant values s = 3 and s = 1.5 
respectively above and below the critical temperature, and a sharp jump at the critical point. 
For finite-size samples the fourth moment is expected to be continuous across the transition, 
but according to the scaling hypothesis the critical value should not depend on the sample 
size if we assume a one-parameter scaling law across the critical point: 

s^s{L/m) (46) 

where L is here the sample length, and ^ is the correlation length which is a function of 
temperature. According to such scaling law s = s{0) at the critical point for any L. 

We have checked this prediction by standard Monte Carlo sampling. For any fixed 
sample size, we have taken a completely random initial configuration, and thermalized it 
at a very high temperature (/? ^ 0.02) by 5 ■ 10*^ complete sweeps. The temperature is 
then decreased by steps of A/? = 0.02. At each step a good thermahzation is achieved 
by 8 • 10^ complete sweeps, and then the averages are evaluated over the successive 2-10^ 
sweeps. Once a sufficiently low temperature is reached (/3 ~ 1), the process is reversed and 
the temperature increased up to the initial value. We have checked that the hysteresis is 
small in all the considered range of temperature. Moreover the small differences observed 
going up and down give a measure of the errors on the configurational averages which have 
been approximated by the mean values. The fourth moment s is reported in Fig.4 for 

19 



= 15, 20, 25, 30. All the curves cross at the same point /3c = 0.53 ± 0.01 as predicted by 
the one-parameter scaling hypothesis. Moreover for very large or very small temperatures 
the correlation length becomes very small and the fourth moment s should approach its 
thermodynamic-limit value s — ^ s(oo) which is expected to be s{oo) = 3 at high temperature 
and s(oo) = 1.5 at low temperature. As shown in Fig. 4 the measured s approaches such 
limits far away from the critical point. 

According to the usual definition of critical exponent 



the scaling equation (^61) allows for an accurate estimate of its value: linearizing s around 
the critical point yields 



where s'(/5c) is the derivative of s as a function of j3. In Fig. 5 a best fit by least squares 
method is reported yielding u = 0.44 ±0.02. Here the error is the statistical one coming out 
from the linear fit. 

Of course this Monte Carlo calculation is far from being the best numerical simulation 
which can be achieved by modern computing machines. Our sample sizes are relatively small 
and a slight shift of the critical point cannot be ruled out. However the estimates for the 
critical temperature and exponent are accurate enough for a comparison with experimental 
findings and for a check of the analytical results of the previous sections, and that is just 
what we needed a the moment. More refined calculations are called for in order to establish 
more accurate predictions. 



Here we summarize and discuss the main findings of the previous sections. According to 
mean-field and finite-size scaling the three-dimensional molecular model has a second order 




(48) 



VI. DISCUSSION 
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continuous transition from an isotropic disordered high-temperature phase to an anisotropic 
two-dimensional ordered low-temperature phase. The d = 3 realization of the model is the 
one which more closely describes real molecular systems. For this reason the d = 3 model 
has been studied by the variational Migdal-Kadanoff method and by numerical Monte Carlo 
simulation. The transition point is characterized by a diverging correlation length according 
to the one-parameter scaling hypothesis which seems to be fulfilled as shown by the data 
of the previous section. On the other hand the d = 2 model is special by itself for its 
equivalence to the two-dimensional Ising model, and for the existence of exact analytical 
results. Thus the d = 3 model can be seen as a non-trivial extension to higher dimension 
of the two-dimensional Ising model. Here "non-trivial" means that the d = 3 molecular 
model does not belong to the universality classes of the standard d = 3 extensions of the 
Ising model (three-dimensional Ising and Potts models). The difference is evident from a 
comparison of the ground state T = configurations: highly degenerate and anisotropic in 
the molecular model (with a two-dimensional character even for higher dimensions); with 
a small degeneration and fully isotropic in the Potts models (including the Ising one as a 
special case). By considering the two-dimensional character of the low-temperature phase, 
the molecular model could be thought to belong to the universality class of the simple two- 
dimensional Ising or three-states Potts models. However in the high temperature unbroken- 
symmetry phase the molecular model is fully isotropic and has a three-dimensional character. 

A formal proof of such statements comes from a comparison of the critical exponents. 
For the three-dimensional molecular model the finite size scaling calculation of the previous 
section yields u = 0.44 to be compared to the two-dimensional two-state (Ising) and three- 
state Potts models whose exponents are u = 1 and u = 0.83 respectively |18|, to the 



three-dimensional Ising model whose exponent is u = 0.64 W^, and to the three-state three 



dimensional Potts model which is known to undergo a first-order transition p3| , p4| . 

The molecular model belongs to a new universality class which is characterized by a sort 
of dimensional transmutation. In fact order takes place in chains which are arranged in 
layers, and the disorder-order transition requires a decrease of the effective dimensionality 
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of the system. In the ordered phase the molecules are correlated inside layers, but there is 
no correlation between molecules which belong to different layers. This understanding of the 
ordered phase is in agreement with our finding that the two-dimensional Migdal-Kadanoff 
variational method for a single layer yields a better prediction for the critical point than the 
three-dimensional method applied to the whole lattice. On the other hand the very same 
two-dimensional variational method provides a convenient analytical tool for describing the 
generic d-dimensional molecular model by a straightforward generalization. 

Prom such arguments the critical point has been given an upper bound by the variational 
method which yields Pc = 0.61, while a lower bound is usually provided by mean- field that 
in this d — 3 case gives Pc — 0.375. The numerical estimate of the previous section Pc — 0.53 
fits nicely inside such bounds. 

Having discussed some formal aspects of the molecular model and few approximate meth- 
ods which throw some light on its phase transition, we would like to make contact with the 
phenomenology. Our main finding - that the order-disorder transition described by the 
model belongs to a new universality class - deserves some experimental test. Transitions of 
this kind have been observed in several systems, as discussed in the introduction. Since the 
critical properties should not depend on the microscopic details of the system we expect that 
the simple molecular model could predict the correct critical exponent of real orientational 
transitions occurring in complex real molecular systems expecially under pressure. New 
experiments are called for in order to test such ideas and explore this broad universality 
class. 
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FIGURES 



FIG. 1. An allowed configuration for d = 2. 

FIG. 2. Ground state configurations for the two-dimensional attractive (a) and repulsive (b) 
models. 

FIG. 3. Numerical solutions of the two-dimensional Migdal-Kadanoff equations for a single layer 
of the three-dimensional molecular model. The critical temperature (3 is reported as a function of 
the field strength h. The stationary point is at /i = /iq = —0.2349 where (3 = (3^. = 0.6122. For 
h > —0.226 there is no physical solution. 

FIG. 4. The fourth moment s = (7^^) versus the inverse temperature (3 for = 15,20,25,30. 
For A = 30 only few points around the critical point have been evaluated. The curves are a linear 
interpolation between points and are reported as a guide for the eye. 

FIG. 5. linear fit for the critical exponent according to Eq.(^). The points have been evaluated 
for A^ = 15,20,25,30. 
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